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Abstract 

The dynamic and static properties of a supercooled (non-entangled) polymer 
melt are investigated via molecular dynamics (MD) simulations. The system 
is confined between two completely smooth and purely repulsive walls. The 
wall-to-wall separation (film thickness), D, is varied from about 3 to about 
14 times the bulk radius of gyration. Despite the geometric confinement, 
the supercooled films exhibit many qualitative features which were also ob- 
served in the bulk and could be analyzed in terms of mode-coupling theory 
(MCT). Examples are the two-step relaxation of the incoherent intermedi- 
ate scattering function, the time-temperature superposition property of the 
late time a-process and the space-time factorization of the scattering func- 
tion on the intermediate time scale of the MCT /^-process. An analysis of 
the temperature dependence of the a-relaxation time suggests that the criti- 
cal temperature, T c , of MCT decreases with D. If the confinement is not too 
strong {D > 10 monomer diameter) the static structure factor of the film coin- 
cides with that of the bulk when compared for the same distance, T — T C (D), 
to the critical temperature. This suggests that T — T C (D) is an important 
temperature scale of our model both in the bulk and in the films. 
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I. INTRODUCTION 



Glass forming materials have been used by mankind since the very early days of civiliza- 
tion. An example are early ceramics from the neolithic period dating back to 5000 B.C. In 
addition to conventional glasses like those used for windows, bottles, etc., polymers represent 
a new class of glassy systems with a large variety of thermal and elastic properties. Due to 
their (in general) low thermal conductivity, polymers are utilized as protective coatings in 
(micro-) electronic devices, optical fibres and other thermally fragile materials |1]-§J. Further- 
more, the high dielectric constant of some polymeric compounds suggests their application 
as porter of electric circuits, a new field under evolution H||. In all these applications, the 
polymer is in contact with another material. This gives rise to an interface. It is therefore 
important to know whether and how the interface alters the properties of the system. In 
particular, one would like to understand to what extent the glass transition temperature is 
influenced by the polymer-substrate interactions. Hence, the study of the glass transition in 
thin polymer films finds a strong motivation from the technological side. 

In addition to the technological importance, the study of the glass transition in thin 
polymer films is also of fundamental interest. It could help to elucidate the nature of 
this phenomenon. According to Adam and Gibbs ||, structural relaxation near the glass 
transition can take place only if many particles move in a correlated way to allow a collective 
rearrangement. The average size of such a "cooperatively rearranging region" defines a 
dynamic correlation length £. The larger £, the smaller the probability of a cooperative 
motion and thus the longer the structural relaxation time are supposed to be. If one assumes 
that £ diverges at some temperature (Vogel-Fulcher-Tammann temperature) below T g , the 
glass transition can be considered as a thermodynamic second order phase transition. 

However, the cooperatively rearranging regions are not directly correlated with static 
density fluctuations. The structure of a glass former changes only slightly upon cooling 
contrary to the hypothesized strong increase of £. As a consequence, scattering experiments 
cannot be used to determine £(T) [0. One has to resort to indirect investigations. A 
possibility consists in studying confined systems. If the system size, L, is finite, one could 
expect that the cooperatively rearranging region cannot grow beyond any bound, its largest 
extension being £ = L < oo. Therefore, the glass transition temperature should depend on 
L and decrease with system size. 

In thin polymer films it is the film thickness, D, which is finite and thus one expects 
T g = T g (D). Indeed, one finds a dependence of T g on D in both experiments fl3Hl3[l and 



computer simulations (14|]. However, the observed change of T g with D strongly depends on 



the system under consideration. If the interaction between the polymers and the substrate is 
attractive, the glass transition temperature of the films becomes larger than the bulk value 
for small film thicknesses ||10|| . Intuitively, this effect can be attributed to chains which are 
close enough to the substrate to 'feel' the attractive interaction. The motion of these chains 
should be slowed down with respect to the bulk. In a thin film almost all chains touch the 
attractive substrate. So, T g should increase. 

On the other hand, measurements (by ellipsometry) of polystyrene (PS) films (of rather 
large molecular weights) on a silicon substrate showed a significant decrease of T g from 
375 K down to 345 K for the smallest film thickness of 10 nm, i.e., a relative change of 10% 
in T g was observed 0. There have also been many experiments in recent years on freely 
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standing polystyrene films (i.e., no solid substrate, but two polymer-air interfaces) [|,|TT|-[T3 
exhibiting a dramatic decrease of T g by up to 20% if the film thickness becomes much smaller 
than the chain size. An interesting explanation of this observation in terms of an interplay 
between polymer-specific properties and free-volume concepts has been proposed |15| . 

This sensitive dependence of T g on the polymer-substrate interaction was also obtained 
from computer simulations [14]. One finds that a strongly attractive wall leads to an increase 
of T g , whereas a weaker attraction has the opposite effect. Therefore, for repulsive walls, a 
decrease of T g can be expected, provided that other parameters of the system, in particular 
the average density, are unaffected by the confinement. As the simulations of our system in 
the bulk were done at constant external pressure [16|;[l7| > we also carried out the simulations 
of the films at constant normal pressure in order to obtain an average density that lies as 
close as possible to the bulk density at the same temperature. This allows to separate the 
effect of confinement from that of the density. 

Note however that the glass transition temperature is an empirical quantity. It is usually 
defined as the temperature at which the viscosity reaches a value of 10 13 poise. So, the 
choice of another number would give a different value of T g . Furthermore, T g depends on 
the cooling rate so that it is not a temperature in a strict thermodynamic sense. Therefore, 
other temperatures have been introduced to characterize the glass transition. One example is 
the so-called Vogel-Fulcher-Tammann (VFT)-temperature, T , at which the system viscosity 
seems to diverge. This quantity is obtained by fitting the viscosity, tj(T), to the empirical 
function, T](T) = i]oexp[E/(T — To)]. An attempt to rationalize the VFT-law is the free 
volume theory |l8|- pl|1 . The main idea of this approach is that a tagged particle can leave 
its initial position only if it finds a "free volume" of size V{>v c in its neighborhood (v c being 
some critical volume of the order of the size of a molecule). It is further supposed that the 
average free volume vanishes at T . Assuming statistical independence of the free volumes 
and using a Taylor expansion of the average free volume around To, one obtains the VFT-law 
for transport coefficients like the viscosity or the (inverse) diffusion coefficient. This brief 
description illustrates that the free volume theory has a phenomenological character. The 
precise meaning of the free volume is unclear and the existence of T , where t>f is supposed 
to vanish, is not derived, but postulated. 

Contrary to the VFT-temperature, the critical temperature of the so-called mode- 
coupling theory (MCT) results from a microscopic approach to the dynamics of supercooled 
(simple) liquids p2|-p4]]. Within the idealized version of the theory, the structural relaxation 
time diverges at a critical temperature T c , while the static properties of the system remain 
liquid-like. This implies that the system vitrifies at T c . Thus, from the point of view of 
MCT, the glass transition is a purely dynamic phenomenon. However, comparisons between 
the theory and experiments p4| , |2"5|| reveal that T c lies in the region of the supercooled liquid, 
where the glass former is only moderately viscous and no absolute freezing occurs. In real- 
ity, there are additional relaxation mechanisms which are not incorporated in the idealized 
MCT and begin to dominate close to and especially below T c . An attempt was made to 
approximately include these relaxation processes in an extended version of MCT P3, but 
the validity of this extension is currently unclear ||26|| . Nevertheless, the idealized MCT 



derives several empirically known phenomena, such as the stretching of the a-relaxation, 
which is often described by the Kohlrausch- Williams- Watts (KWW)-law (oc exp[— (t/r) 13 ], 
where /3<1), or the time-temperature superposition principle. Furthermore, it also makes 
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new predictions, such as the space-time factorization property, which have been tested both 



in experiments |24|,|25| and computer simulations 

An application of MCT to our polymer model in the bulk |IT| , |28| -|3Tf revealed that the 
theory represents a suitable framework to analyze the dynamics of the supercooled melt. 
Therefore, we attempt to test whether MCT can also be applied to the supercooled polymer 
films. We find that even for an extremely thin film of three monomer layers, some features of 
the dynamics at low temperatures can be successfully described by the MCT. This analysis 
yields T c as a function of film thickness and shows that T c decreases with D. Within the 
error bars the Vogel-Fulcher-Tammann temperature exhibits the same D-dependence so that 
we also expect the glass transition temperature of our model to behave analogously due to 
Tq(D) < T g (D) < T C (D). 

The paper is organized as follows: After a presentation of the model in section [Tl], sec- 
tion |H focuses on the influence of the confinement on static properties of the system. A 
discussion of the dynamics is the subject of section |TV| . In section |V| we determine the 
dependence of T c and T on film thickness and the last section summarizes our conclusions. 



II. MODEL 



We study a Lennard- Jones (LJ) model for a dense polymer melt []16],|32] of short chains 
(each consisting of 10 monomers) embedded between two completely smooth, impenetrable 
walls 



3.34 



Two potentials are used for the interaction between particles. The first one is 
a truncated and shifted LJ-potential which acts between all pairs of particles regardless of 
whether they are connected or not, 



LJ-ts \ 







U LJ (r c 



if r < r c , 
otherwise , 



where U^j(r) =4e (cr/r) 12 — (cr/r) 6 and r c = 2 x 2 1 / 6 a. The connectivity between adjacent 
monomers of a chain is ensured by a FENE-potential [ |32|| . 



k 



1 - 



Rc 



(1) 



where k = 30e/<7 2 is the strength factor and Rq = 1.5<r the maximum allowed length of a 
bond. The wall potential was chosen as 



(2) 



where z = |-z pa rticie — Zwa\\\ (zw&a = ±D/2). This corresponds to an infinitely thick wall 
made of infinitely small particles which interact with inner particles via the potential 
45e(o"/r) 12 /(47rp wall (T 3 ) where p wa n denotes the density of wall particles. The sum over the 
wall particles then yields e(a/z) 9 . All simulation results are given in Lennard- Jones (LJ) 
units. All lengths and energies are measured respectively in units of a and e, temperature 
in units of e/k-Q (/cb = 1) and time in units of (mcr 2 /e) 1 / 2 . 
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The left panel of Fig. |T| compares the bond potential, i.e. the sum of LJ- and FENE- 
potentials, with the LJ-potential. It shows that the bonded monomers prefer shorter dis- 
tances than the non-bonded ones. Thus, our model contains two intrinsic length scales (see 
right panel of Fig. [I] for a schematic illustration). Since these length scales are chosen to be 
incompatible with a (fee or bec) crystalline structure and since our chains are flexible (no 
bond angle or torsion potentials), one could expect that the system does not crystallize at 
low temperatures, but remains amorphous j3~5 . 



This expectation is well borne out, as Fig. |2| illustrates. The upper part of the figure shows 
a snapshot of a film of thickness D = 20 and at temperature T = 0.44 which corresponds to 
the supercooled state. The visual inspection of this configuration suggests that the structure 
is disordered. This is corroborated by an analysis of the static structure factor S(q). 

The lower part of Fig. |2| shows S(q) for a film of thickness D = 10 at temperatures 
corresponding to the normal liquid state and to the supercooled state. S(q) is calculated 
parallel to the walls [i.e. q— \q\, q = (q x , q y )] by averaging over all monomers in the system. 
Figure |2] shows that the maximum value of the S(q) increases at lower temperature and 
that its position is slightly shifted towards larger g-values. Thus, the average interparticle 
distance decreases with decreasing temperature, since the density of the film increases in our 
simulations at constant pressure However, these structural changes are rather small. 
Even at a very low temperature of T = 0.42 which is quite close to the critical temperature 
of the system at this film thickness (T C (D = 10) ~ 0.39), no qualitative difference is observed 
between the structure factors at low and high temperatures. The packing of the system thus 
remains liquid-like (i.e. amorphous) at all studied temperatures. 

All simulations have been carried out at constant normal pressure Pn ext = P = 1 p4 



However, to adjust the normal pressure, we do not vary the wall-to-wall separation, D, but 
the surface area. For each temperature, the average surface area is calculated by an iterative 
approach |j36| . The system is then propagated until the instantaneous surface area reaches 
the computed average value. At this point the surface area (and thus the volume) is fixed 
and a production run is started in Al/T-ensemble, where the temperature is adjusted using 
the Nose-Hoover thermostat [37. 3^] . This thermostat slows down or accelerates all particles 



depending on the sign of the difference between the instantaneous kinetic energy of the 
system and the desired value imposed, i.e. 3iV/c B T/2 (N is the number of particles) |37|-442 



One may therefore ask how reliable the resulting dynamics is when compared to pure 
Newtonian dynamics in the microcanonical (NVE) ensemble. This question was already 



examined for the bulk system in Ref. and for our confined model in Refs. P3"| , fP| . In 
both cases, the results obtained from constant energy (NVE) simulations and from the 
Nose- Hoover thermostat are identical. More details about the applied simulation techniques 
can be found in p4],|36|,^4 1 . 



III. STATIC PROPERTIES 

In this section, we want to discuss the influence of confinement on both the chain confor- 
mation and the static properties of the melt. When computed within thin layers parallel to 
the wall, the structure of a chain varies with the distance from the wall. On the other hand, 
the chain structure averaged over the whole film does not depend much on film thickness. 
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Contrary to this insensibility of the average chain conformation to the confinement, the 
dense packing of the melt exhibits a pronounced dependence on D. 

A. Effects of the Confinement on the Structure of a Chain 

Let us first look at the T-dependence of the single chain structure factor S c (q). Simi- 
larly to S(q), S c (q) is also calculated for q- vectors in direction parallel to the walls and by 
averaging over all chains in the system. 

The upper panel of Fig. |3] shows S c (q) for a film of thickness D = 20 at two representative 
temperatures: a relative high temperature corresponding to the normal liquid state (T = 
1) and a low temperature representative of the supercooled state (T = 0.44). Contrary 
to the structure factor of the melt (see Fig. 0), S c (q) is practically independent of the 
temperature, not only on scales larger than the radius of gyration, but also on the local 
scale of the intermonomer distance. This may be rationalized as follows: The interactions of 
the monomers along the backbone of a chain do not include specific potentials for the bond 
or torsional angles, which could make the chain expand with decreasing temperature. This 
would lead to a much stronger temperature dependence than that resulting from potentials 
used. Since the bond potential is very steep around the minimum and is in general very 
large, the bond length is essentially independent of T (b = y/ (b ■ b) = 0.966 at T = 1, 
b = a/ (b ■ b) = 0.961 at T = 0.46. Here, b is the bond vector). The main effect is that 
the overall size of a chain slightly shrinks with decreasing temperature because the density 
increases. This also leads to a weak increase of the peak of S c (q) at q = 7.6, which is, 
however, not visible on the scale of Fig. |3] and negligible compared to that of S(q) in Fig. |2|. 

Furthermore, Fig. |3| shows that, for q < 27r/i?g Ulk , the g-dependence of S c (q) can be 
described well by a Debye-function, 

S^'il) = ^pgy («P(-^ D *°) " 1 + 1 2 "l D,b,e ) • (3) 

with Debyc _ ^>2 buik^ wnere ^2 bulk wag ca j cu i a ted independently and inserted in Eq. (|3]) 
[N p is the number of the monomers of a chain (degree of polymerization)]. The same 
observation has already been made for the model in the bulk |TjJ. This agreement may be 
expected because the scattering for q < 27c/R^ nlk is determined by the overall size of a chain. 
On this scale, the difference between the freely-jointed chain model utilized to derive the 
Debye function and the actual model of the simulation does not matter. 

The lower panel of Fig. |3] shows S c (q) for films of various thicknesses D at T = 1. 
Obviously, the chain conformation averaged over the whole film does not depend much on 
D. Contrary to the small effect of temperature on S c (q), the independence of film thickness 
is less intuitive. Figure [| illustrates this point. In the upper panel we compare the parallel 
and perpendicular components of the radius of gyration and of the end-to-end distance, 
i?g || , Rg t ± and R^ e m, R^^, for D = 20 and T = 1 (normal liquid phase). The components 
are plotted versus the distance, z, from the (left) wall, where z denotes the position of 
chain's center of mass. So, R 2 g 11(2), for instance, is the radius of gyration parallel to the wall, 
which is averaged over all chains whose centers of mass are located at z. The figure shows 
that both the radius of gyration and the end-to-end distance agree with the bulk value if 
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z > z w + 2_Rg Ulk . Here, z w « 1 denotes the physical position of the wall, i.e., the smallest 
distance between a monomer and z wa \\ = D/2. As the chain's center of mass approaches 
the wall, R? N and R 2 n first develop a shallow minimum and then increase to about twice 
the bulk value followed by a sharp decrease to zero in the very vicinity of the wall where 
practically no chain is present. On the other hand, the perpendicular components, Rg± 
and R 2 e ±, fi rs t pass through a maximum before decreasing to almost at the wall. This 
behavior has been observed in several other simulations (see EH] and references therein), 
also for larger chain length than that used here |46| , |47[| . 

These results suggest that a spatially resolved version of the chain structure factor, 
S c (q, z), should depend upon the position of the chain with respect to the wall. To test this 
idea we divide the system into layers of thickness 1/4 (in units of a which roughly corresponds 
to the monomer diameter) and evaluate S c (q, z) by taking into account all monomer pairs 
within the layer which belong to the same chain. Here, z is the distance of the middle of such 
a layer from the wall. The lower panel of Fig. |3] shows S c (q, z) for an extremely thin film of 
thickness D = 5 at T = 0.35 and for two typical layers: a layer in the film center (z = 2.375) 
and a layer close to the wall (z = 1.125) [note that all layers with a distance smaller than 
z = 1.125 to the wall are practically empty. Note also that there is no layer whose middle 
lies exactly in the film center z = 2.5. Only the boundary of the central layer "touches" 
the film center so that its middle is closer to the wall than D/2 = 2.5]. In the center of the 
film S c (q,z) can be described by a Debye function with R 2 Deb y e = 1.93, whereas a larger 
radius (R 2 Debye = 2.4) must be chosen close to the walls. These values for R 2 Deb y e are taken 
from the profile of the radius of gyration for D = 5 at T = 0.35, which yields R 2 = 2.4 and 
Rg = 1-93 when averaging over the intervals 1 < z < 1.25 and 2.25 < z < 2.5, respectively. 
Furthermore, Fig. |] shows that S c (q,z) for the film center coincides with S c (q) obtained 
by averaging over the whole system. Thus, the local variations of S c (q, z) close to the wall 
disappear when the whole system is considered. Qualitatively, this point can be understood 
by noting that there are only few chains close to the wall, where R 2 Jz) is larger than the 
bulk value, while most of chains are in the inner part of the film, where R 2 m (z) is close to 
R 2 bulk . Averaging over the whole system thus tends to cancel the effects of the walls on the 
chain conformation. 

This point is further examined in Fig. [|. The upper panel of the figure shows 1-5-Rg n, 
3i?g j_, averaged over the whole film, and R 2 = R 2 i, + R 2 ± versus temperature for two film 
thicknesses D = 5 and D = 20. For all temperatures the parallel component is larger than 
the bulk value, whereas the perpendicular one is smaller. The disparity between the two 
components becomes more pronounced for small film thickness. Nevertheless, the sum of 
these quantities, R 2 ., is fairly close to the bulk value for both film thicknesses. It slightly 
shrinks with decreasing temperature, since the density of the films increases. Similar trends 
are observed if D is varied and T is kept constant. The lower panel of Fig. |5| depicts 1-5-Rg n, 
3i?g ± and R 2 as a function of film thickness at T = 0.5. R 2 depends only weakly on film 
thickness and is close to the bulk value. This explains why the chain structure factors of 
Fig. H are essentially independent of D. 
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B. The Packing Structure of the Melt 



Figure ^| compares the structure factor of the melt, S(q), in the bulk with that of films 
of various thicknesses at T = 1 (normal liquid state; upper panel of the figure) and at 
T = 0.46 (supercooled state; lower panel of the figure). Qualitatively, the behavior of the 
films and the bulk is identical. The structure factor is small at small g-values, reflecting 
the low compressibility of the system. Then, it increases and develops a peak at g max which 
corresponds to the local packing of monomers (27r/g max ~ 1) before it decreases again and 
begins oscillating around 1, the large-g limit of S(q). This behavior is characteristic of dense 
amorphous packing. 

However, there are quantitative differences which become more pronounced at low tem- 
perature: While S(q) of the film of thickness D = 20 (almost) coincides with the bulk data 
at T — l, deviations are clearly visible at T = 46. Quite generally, the most prominent differ- 
ences between the bulk and the film are found for small q and for g max . The compressibility 
of the film is higher, the value of g max is shifted to slightly lower q and the magnitude of 
S(q max ) is smaller than in the bulk. Keeping the film thickness fixed, one can observe similar 
changes of S(q) as the temperature increases (see Fig. |2|). Therefore, the local packing of 
the monomers in the films seems to resemble that of the bulk at some higher temperature. 
Since the local structure of the melt has an important influence on its dynamic behavior in 



the supercooled state p2| , Fig. g suggests that the film relaxes more easily than the bulk at 
the same temperature. Indeed, we will see later that the dynamics of the system is much 
faster in the film than in the bulk when compared at the same temperature. 



IV. DYNAMICS 

This section discusses the dynamics of the films at low temperatures and for various 
thicknesses, D, ranging from about 3 to about 14 times the bulk radius of gyration. To this 
end, the incoherent intermediate scattering function and various mean-square displacements 
were calculated. We will see that, despite the geometric confinement, the films exhibit 
several dynamic features which are in agreement with predictions of mode-coupling theory 
(MCT) [||^5|. In this respect, the films behave as the bulk |I7|,g?HDl. However, the onset 
of MCT effects is shifted to lower temperatures compared to the bulk. The presence of the 
smooth walls accelerates the dynamics and this influence of confinement is the stronger, the 
smaller D. 



A. Confinement leads to Faster Dynamics 

An interesting dynamic correlation function is the incoherent intermediate scattering 
function (f) q (t). It measures density fluctuations on various length scales, which are caused 
by the displacement of individual particles. For a planar system, we define <jf q (t) for q- vectors 
parallel to the wall, i.e., 

1 N 

W = (ivE ex P fall Kll(*) -r<.||(0)])) • (4) 
i=i 



S 



Here, iV is the total number of monomers in the system, q« = (q x ,q y ), q= \qu \ = y/q% + q'y 
and r it \\ = (xi,yi). 

Figure [7] compares the relaxation of <fi s q (t) at the maximum of the static structure factor 
(g = 6.9) in the bulk with that in films of various thicknesses ranging from D = 5~3.5R g to 
£> = 20«14R g . The temperature studied, T = 0.46, is slightly above T c bulk (=0.45) 



In the bulk, a two-step relaxation is observed: At very short times, (jf q (t) can be described 
by (ffJi) = 1 — (Q s t) 2 /2 with the microscopic frequency Q s = qy/k^T [52]. This corresponds 



q 

to free particle motion. At later times, the relaxation of (jf q (t) is strongly protracted. There 
is an intermediate time window (/^-relaxation regime of MCT ||22||), where the correlator 
changes rather slowly with time before the final structural relaxation (a-relaxation) sets 
in at long times. This two-step decay is a characteristic feature of 4> s q (t) for temperatures 
close to T c and reflects the temporary "arrest" of a monomer in its local environment ("cage 
effect" It has been analyzed in detail in Refs. |§|13,|8|,|9|. 



Compared to the bulk, the film data relax faster. If the film thickness decreases, this 
acceleration is enhanced and the two-step decay gradually disappears. At D = 5, no interme- 
diate /3-relaxation is observed. The same changes also occur in the bulk if the temperature 
increases |IT| , |2"P| , f4"g| . This suggests that the inverse film thickness qualitatively plays a similar 



role as the temperature. 

The acceleration of the dynamics in the film compared to the bulk is not limited to the 
main peak of S(q). It is also found for other q- values, for instance, for very small wave 
vectors. Since the time dependence of <fi B q (t) is directly related to the monomer mean-square 
displacement (MSD) in the low-g limit, it is instructive to illustrate the acceleration of the 
dynamics by an investigation of the MSD. For a polymer system, various kinds of MSD's 
may be defined. An important example is the mean-square displacement of the innermost 
monomer, 

M 

^w = ^E(i r rw- r ™ r (°)] 2 ) • ( 5 ) 
i=i 

Here, M is the total number of chains in the melt and r^ er = (x* nner , yl nner ) is the projection 
of the position vector of the innermost monomer of the i-th chain onto the xy-plane (which 
is parallel to the walls). The factor 3/2 is introduced to simplify the comparison with bulk 
results. It takes into account that only two independent directions (x and y) contribute to 
the MSD in the film, whereas all three directions are considered in the bulk[} Similarly, one 
defines the mean-square displacements of the chain's center of mass, 

fliW = 2M E (TO*) " ^ ' ^ 
where is the projection of vector to the i-th chain's center of mass onto a plane parallel 



As most of the comparisons with bulk results are done using the film data in parallel direction, we 
drop the index "||" to simplify the notation. To avoid ambiguities, when film data in perpendicular 
direction are discussed, an index "_L" will then be used. 
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to the wall. Again, the factor 3/2 accounts for the difference in the number of independent 
components, as above. 

Figure [8] depicts gi(t) for two typical temperatures: T = 1 (upper panel), which cor- 
responds to the normal liquid state, and T = 0.46 (lower panel), which belongs to the 
supercooled state in the bulk. Both panels compare gi(t) for the bulk and for films of var- 
ious thicknesses. The influence of the walls is rather small at T = 1 so that gi(t) of the 
bulk almost overlaps with that of the film if D > 10. However, the lower panel of Fig. |§| 
shows that the effect of the walls on the mobility becomes significant at all studied thick- 
nesses with progressive supercooling. Outside the initial ballistic regime (gi(t) = 3Tt 2 ), the 
motion resembles that of the bulk, but is the faster, the smaller the film thickness. For the 
bulk and D > 7, g%(t) exhibits several regimes. In agreement with the predictions of the 
mode-coupling theory fl22|-f2"4[| , a plateau regime emerges after the ballistic motion. At low 
temperature, the tagged particle remains temporarily in the "cage" formed by its neighbors. 
However, contrary to simple (atomic) liquids, where a direct crossover from the plateau to 
the diffusive regime occurs, an intermediate subdiffusive regime emerges due to the connec- 



tivity of the monomers ||50|| . In this regime, which is present for all D, the center of mass 
already crosses over to the asymptotic diffusive motion, g%(t) — t, whereas the motion of 
the innermost monomer is described by a power law gi(t) ~ t x with an effective exponent 
x ~ 0.63. The innermost monomer reaches the center of mass only if g\(t) is larger than 
the end-to-end distance of a chain. In this limit, the diffusive motion is dominated by the 
motion of the chain's center of mass and gi(t) coincides with gz{t). 

In addition to the parallel displacements the analysis of the motion perpendicular to 
the wall is also interesting because the oscillations of the monomer density profile, which 
can propagate through the whole film for small D and low T, could possibly suppress 
perpendicular motion substantially. To check that, Fig. |9] shows the MSD of all monomers, 
go(t), computed in direction parallel to wall, 

3 N 

= ^£<Ku(*)-rd 2 ) , (7) 

and perpendicular to it 



2N 



3 - 

Here, N is again the total number of monomers and n = {xi, yi). Furthermore, Zi denotes 
the z position of z-th monomer. The factors 3/2 for go(t) and 3 for <7o,±(*) account for 
the difference of independent components: two in parallel and one perpendicular direction 
compared to 3 in a bulk system. 

The upper panel of Fig. |] depicts go(t) at T = 0.46 (supercooled state). It compares 
the bulk data with the displacements parallel and perpendicular to the walls in films of 
thicknesses D = 5, 7 and 20. The comparison of go,±{t) for the films and of go(t) for the bulk 
reveals that the confinement does not only accelerate the dynamics in parallel, but also in 
perpendicular direction if go,±(t, D) is sufficiently smaller than the film thickness. For all D, 
one then finds go t ±(t,D) > #o ulk (£) and go t ±(t, D{) > go,±(t,D2) if D\ < D 2 . However, this 
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increase of the mobility is less pronounced than that of the dynamics in parallel direction, 
so that, for a given film thickness, g 0t \\(t,D) > g 0> ±(t,D). 

Of course, the perpendicular displacement cannot grow infinitely. It must be limited by 
the film thickness. This is illustrated by go 7 ±(t,D = 5) which crosses over to a constant of 
approximately 5 at late times. In fact, using the density profile, p(z), one can compute the 
large time limit of go,±{t) by 

-D/2 r+D/2 

dz / dz'p{z)p{z'){z-z'f 
= 3 ^ /2 J -f +D/2 • (9) 

K J-D/2 J 

Equation (Q) can be understood as follows: The contribution of two points z and z' to the 
long time limit of go,±(t) is equal to [z— z') 2 multiplied by the probability p {2 \z,t] z',t') dzdz' 
of finding a tagged particle at a given time, t', in the interval [z', z'+dz'] provided that it was 
at a previous time, t in [z, z + dz]. As \t — f\ grows, this probability becomes independent of 
the initial position of the particle, i.e. p( 2 \z,t; z' ,t')=p(z,t) p{z' ,t'), where p(z, t)dz is the 
probability of finding a particle at time t in [z, z + dz]. Obviously, this probability does not 
depend on the choice of the time origin, i.e. p(z, t) —p(z) —p(z)/ f p(z')dz'. Putting all this 
together and adding a factor of 3 for the sake of comparison with bulk results, we obtain 
3(z — z') 2 p[z) p(z') dzdz' j (J p(z)dz) 2 for the contribution of the pair (z, z') to linu^oo go,±(t). 
Integration over both variables z and z' yields Eq. (^). 

For a film of thickness D = 5 at a temperature of T = 0.46 insertion of the density 
profile (obtained from the simulation) in Eq. (H) yields lirn^oo go,±{t) =4.667. This is the 
value to which go,±(t,D = 5) converges for large t (Fig. |9p. Note that by replacing the 
monomer density profile p(z) in Eq. ([|) by the density profile of the innermost monomer, 
one can obtain the long time limit of gi,±{t), i.e. 5 , i i j_(t = oo). In a similar way, Eq. (Q) can 
be adapted to compute g^± (t = oo) by using the density profile of end monomers and/or 
(73 _l (t = oo) by replacing p(z) by the profile of the density of the chain's center of mass, 

Pern (z) ■ 

The general validity of Eq. ([§) is demonstrated in the lower panel of Fig. ^. Here, 
9o,±(t) an d gi,±(t) are plotted for two temperatures, T = 0.46 and T—l. The corresponding 
long time limits have been computed as discussed above. At both temperatures the results 
obtained by Eq. @ are in good agreement with long time behavior of go,±(t) and gi,±(t). 
The inset of the lower panel of Fig. ^ depicts the density profiles of all monomers and of the 
innermost monomer at both investigated temperatures. As seen from this inset, the density 
close to the walls increases appreciably at lower temperatures. Furthermore, the formation 
of strong density peaks close to the walls is accompanied by a slight decrease of the density in 
the center of the film. As a consequence, the relative weight of larger transversal distances in 
Eq. ([]) increases compared to that of smaller [z — z'). One can therefore expect an increase 
of 9o,±{t = oo) at lower T. As the density of the innermost monomer changes in a similar 
way with temperature, we also expect an increase of gi t ±{t = oo) when decreasing T. This 
expectation is nicely born out in the lower panel of Fig. Rl 
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B. Incoherent intermediate scattering function 



The discussion of the previous sections suggested that mode-coupling theory (MCT) 
could also be a relevant theoretical framework to describe the dynamics of the supercooled 
polymer films. In this section we want to test this suggestion by an analysis of the incoherent 
intermediate scattering function (jfJt). 

A quantitative application of MCT to the simulation data requires an intricate fit proce- 
dure which must simultaneously optimize several parameters subject to various theoretical 
constraints (some of them have to be independent of temperature, others independent of the 
wave vector). Before attempting this analysis simple tests should be carried out to check 
whether the approach is worthwhile at all. Two such tests can be performed. 

Mode-coupling theory predicts that there is an intermediate time window in which the 
scattering function remains close the time- and temperature independent non-ergodicity 
parameter f^ c . This time window is called /3-relaxation regime |2"2"| , |2"5| . In this regime the 
scattering function can be written as 

m = it + KGit) , (io) 



where G(t) and h s q are the /3-correlator and the critical amplitude, respectively P2] , |23 
Equation ([K]) shows that the time-dependent corrections to /| c have an important property. 
The space- and the time dependences factorize from one another. 

This "factorization theorem" [0,^] suggests a simple test p^f49] , |^|52"|| which uses the 



simulation data directly without invoking any fit procedure. If £' and £" denote two times 
belonging to the /3-regime, then the ratio 



<g) ~ ggO Git) - G{t') 
cjfJt") - Mt>) G(t") - G(f) 



RS(t) = ^ \ = „>" } -y \ = R(t) (11) 

<? v ' A&(+"\ -LsU/\ ni4.n\ rifj-i\ w 



should only depend on temperature and time, but not on q, provided Eq. ([10]) holds. 

Note that (j) s „(t) varies slowly for times around the plateau. Therefore, the denominator 
of Eq. ([□]) is fairly small and the accuracy of the test is predicated upon an appropriate 
choice of the parameters £" and if. To obtain a satisfactory signal-to-noise ratio, one would 
like to make t" and tf as different as possible. However, one has to be careful not to take t" 
and £' outside the /3-relaxation regime, where Eq. (fLl|) is no longer valid. We find that t" = 1 
and £' = 50 is a reasonable compromise for all studied film thicknesses D = 5, 10 and D = 20. 

With this choice two observations can be made from Fig. |5|: First, there is indeed an 
intermediate time window for all D where the correlators, measured at different q, collapse, 
whereas they splay out at both shorter and longer times. This is a qualitative evidence for the 
factorization theorem. However, the best agreement with Eq. (|TT1) is obtained for D— 10. In 
the other two cases the superposition of the scattering functions for £" < £ < £' is not as good. 
A possible reason for this difference could be that, for D— 10, Eq. fllTD is tested at a distance 
of T — T C (D = 10) = 0.42 — 0.39 = 0.03 from the critical temperature of this film thickness. 
The tests for _D = 20 and D = 5, however, correspond to T — T C (D = 20) = 0.46 - 0.415 = 0.045 
and T — T C {D = 5) = 0.35 — 0.305 = 0.045, respectively. As Eq. (|TT| ) is an asymtotic relation 
which is expected to hold the better, the closer T is to T c , Fig. ^| still suggests that the 
factorization theorem is not only satisfied in the bulk [0,H3], but also in the polymer films. 
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The second observation concerns the order of the g-values before and after the /3-regime. 
This order is preserved. The top curve, (jf ^t), at short times is also the top curve in the 
a-regime. Similarly, the bottom curve, 0^=18.5^), before the /3-regime also remains below all 
other g-values when leaving the /3-regime again. This behavior reflects the theoretical pre- 
diction that the short- and long-time corrections to Eq. ( |i~0"D exhibit the same g-dependence. 



It was pointed out in simulations of a binary Lennard- Jones mixture and also found for 
our model in the bulk | |48| | . 

The second test of the applicability of MCT deals with the late-time relaxation of <fi s q (t). 
An important prediction of the theory for the a-process is the time-temperature super- 
position principle (TTSP) P2|p51 - This means that <p q (t) is not a function of time and 
temperature separately, but only of the scaled time t/r q (T), where r q (T) is the a-relaxation 
time. So, we have 

<P q (t) = F q (t/r q ) . (12) 

The function F q (t/T q ) can often be well approximated by a Kohlrausch- Williams- Watts 
(KWW) function 

F q {t/rf) = ffe W [{t/rfff] . (13) 

For time-temperature superposition to hold the amplitude and the stretching exponent 
0q- must be independent of temperature. Only the Kohlrausch relaxation time is a 
function of T. MCT predicts that is proportional to the a-relaxation time and increases 
upon cooling as p2] , |25[| 

T?KT q ~(T-T c )^ . (14) 



Equation ([14] ) also implies that any time from the window of the a-process should exhibit 
the same temperature dependence and can thus be used to test the TTSP. This means that 
it should be possible to collapse the late time decay of the scattering functions, measured 
at different T, onto a common T-independent master curve by plotting <fi q {t) versus tjr q . A 
convenient definition of r q is to simply read off the time when <f) s q (t) has decay to a certain 
value. Again, such a test has the advantage that no complicated fit procedure is involved. It 
works directly with the simulation data. A possible choice is 4> s q {T q ) =0.1 This low value 



warrants that the scattering function has decayed sufficiently so that possible perturbations 
from the /3-relaxation are completely negligible. 

The resulting master curves are shown in Fig. [H] for D = 5, 10, 20. For all film thicknesses, 
even for the extreme case of D = 5, which corresponds to three atomic layers only, the TTSP 
is borne out by the simulation data and extends to shorter rescaled times with decreasing 
temperature, as predicted by MCT for homogeneous systems. In this respect, the films 
behave identical to the bulk fi~7| , f4"9"|| . However, the results for D = 10 suggest that there 



could also be a qualitative difference. They include T = 0.4, which is very close to the 
estimated T c , i.e. T — T c = 0.01. Contrary to the bulk [ |i^ , [49|1 , the film data at this T — T c 
exhibit no apparent violation of the TTSP. Whether this is a general property of the confined 
systems or just a special feature of D = 10 is not clear at present. 

While the scattering functions exhibit time-temperature superposition for all thicknesses 
studied, they do not superimpose if different thicknesses are compared. This is illustrated in 
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Fig. |12[ The figure shows (jf versus t/r q for D = 5, 20 and the bulk at comparable distances 
to the corresponding critical temperature (i.e., T — T c ~ 0.045). Obviously, the shape of 
the late-time relaxation depends on D. Note that, compared to j3^(D = 20), the stretching 
exponent of the thinner film, i.e., 0f{D = 5) is closer to that of the bulk. A possible 
explanation could be as follows: Due to the presence of the walls, there is a distribution of 
the relaxation times along the transversal direction. Regions closer to the wall decay faster 
and thus exhibit a smaller relaxation time. On the other hand, as the temperature decreases, 
the influence of the walls is "felt" throughout the film for all thicknesses studied. There is 
practically no region of (bulk-like) constant relaxation time in the inner part of the film. 
Now, if we divide a film into layers in which the relaxation time is calculated, the number 
of layers is smaller in a thin than in a thick film. The simplest assumption then is that the 
larger number of layers in the thick film leads to a broader distribution of relaxation times. 
Of course, if D increases further, there will be a region of bulk-like behavior which grows 
in the middle of the film and starts to dominate properties of the system for very large D. 
The film thicknesses studied here, however, are far from this limit. 



V. DETERMINATION OF T C (D) 



The analysis of the preceeding section showed that the relaxation of the films is qual- 
itatively compatible with predictions of mode-coupling theory and that it speeds up with 
increasing confinement. To some extent, the dynamics of the films corresponds to the be- 
havior of the bulk obtained at a higher temperature. This suggests that confinement reduces 
the characteristic temperatures of the film compared to those of the bulk. In this section, 
we want to quantify this reduction by determining the critical temperature T C (D) and the 
Vogel-Fulcher-Tammann temperature T (D). 

711 predicts that the relaxation times of any correlation func- 
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tion, which couples to density fluctuations, should exhibit the power-law temperature de- 
pendence of Eq. ([14]) if T is close to T c . Indeed, the a-relaxation times of the incoherent and 



coherent scattering functions in the bulk may be described by Eq. (|T4[) in some g-dependent 
temperature interval fl?,4"9|| . Deviations are found both for temperatures too close and to far 
away from T c bulk . These deviations are theoretically expected. Since Eq. (|14]) is an asymtotic 
result of the idealized MCT, it can be violated if T is very close to T c due to relaxation chan- 
nels which are not treated by the idealized theory, and very far from T c where the asymtotic 
formula is no longer applicable. Similarly, the g-dependence of the temperature interval, 



where Eq. fll4) is valid, has been rationalized by calculating the leading-order corrections to 
the asymptotic behavior within the framework of idealized MCT [ 5B| , |51 1 . Due to our findings 
in the bulk and due to the results of the previous section we expect to obtain similar results 
when analyzing the films. 

To examine this expectation, we define relaxation times as the time needed by a given 
mean-square displacement, like g\, g 3 or g [see Eqs. @, (§) and (0)], to reach the monomer 
size 

gi(t = r) := 1 (defining equation for r) . (15) 
This is a reasonable choice because the bulk analysis revealed that a monomer, which has 



moved across its own diameter, contributes to the a-relaxation 48 . Using Eq. (jl5|) we 
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computed r{gi = 1) as a function of temperature for various film thicknesses, where, in 
addition to g$, g\ and g%, the MSD of the end monomers, 



M 



(16) 



has also been used. The resulting relaxation times were then fitted by Eq. (|14]) as follows: 
First, all (three) parameters of Eq. were adjusted. The results for 7 obtained from 
r(g = l), r(gi = l), r(g 3 = l) and r(g 4 = l) agreed well within the error bars. Therefore, we 
fixed 7 at the average value for the given film thickness and repeated the fits to optimize 
the remaining parameters. 

Table | compiles T C (D) and 7msd(-D) obtained in this way. Figure [13] shows a representa- 
tive example of this analysis. The upper panel depicts r -1 / 7 versus T for a film of thickness 
D = 5. The intersection of the straight lines (MCT-fit results) with the T-axis yields the 
critical temperature at this film thickness: T C (D = 5) = 0.305 ± 0.005. Note that, despite 
the highly non-linear relationship between the MSD's used to define the various r's, all fits 
yield the same T c . 

To test this analysis the determined critical temperature can be used to linearize the 
relaxation time by plotting r versus T — T c on a log-log scale. The lower panel of Fig. O 
shows that the power law ([14]) is a good approximation close to T c . Analogous to the bulk, 



there are deviations for large T— T c , where the asymptotic regime is left. On the other hand, 



no deviations are observed for small T — T c . In the bulk [|TT[p9"| , they were only found for 
T—T c < 0.02 and thus for temperatures smaller than those simulated for D = 5 up to now. 

Furthermore, as indicated by the solid line in the lower panel, r(T) can also be described 
by a Vogel-Fulcher-Tammann (VFT) equation, i.e., by 



t(T) oc exp 



c{D) 



T-T (D) 



(17) 



where c is a constant which can depend on film thickness. The possibility of describing the 
same data by both a power law (MCT) and a VFT-fit has also been observed for the bulk 
(see Fig. 10 in |16]). We therefore use the VFT-formula as an independent approach to 



determine the variation of T with film thickness. Table | contains the results. A plot of 
T C (D) and T (D) is shown in Fig. [14]. Since we expect T < T g < T c , the figure suggests 
that also the glass transition temperature, T g , should be reduced for stronger confinement. 
Qualitatively, this result parallels those reported for experiments on supported || and on 
freely standing polystyrene films as well as for MD simulations of a polymer model 

similar to ours in the case of weak monomer-substrate attraction H . 

The critical temperatures, T C (D), were determined, for instance, from the mean-square 
displacement of all monomers [see Eq. ([IB])] and thus from a quantitiy which corresponds to 
the low-g limit of the incoherent scattering function. Figure |15| shows that the same T C (D) 
can also be used to linearize the relaxation times r q at maximum of the static structure 
factor. Here, r q was defined by 4> s q (t = r q ) = 0.3. The upper panel depicts a log-log plot of 
r q versus T — T c for D = 5,7, 10, 15, 20 and for the bulk. It illustrates that a shift of the 
temperature scale by T C (D) leads to a very good superposition of the film data for all T 
studied and also of the film and the bulk if T — T c > 0.1. For smaller T — T c deviations occur. 
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The relaxation time of the bulk increases less steeply than that of the films. This implies 
that the exponent 70 of Eq. (14) is smaller in the bulk than in the films. If one fits Eq. (14) to 
that T-interval, where the film data are linear, one obtains 7^-values which range between 
70 (-D = 10) = 2.68 and ^(D = 5) = 3.24 and which do not increase monotonically with 
decreasing D [see table M. The lower panel of Fig. [15| illustrates the same data in a slightly 
different way. Here, r q 1 ^ 14 ' is depicted versus T — T C (D). The exponents, 7^, used for this 
plot are obtained from fits to Eq. ( |I4]) using r q as input data and keeping T C (D) constant. As 
expected, curves at all film thicknesses converge towards the origin of coordinates. In other 
words, the computed critical temperatures are consistent with a temperature dependence 
of Tg according to Eq. ([14]). On the other hand, the curves splay out at large T — T C (D) 

indicating that the prefactor, B, in r q 1 =B \T — T C (D)\ depends on the film thickness, D. 

As mentioned above, although not monotonically, the exponent 7^ seems to increase 
with decreasing D [see also table 0. To some extent, this finding is unexpected, given the 
variation of the Kohlrausch exponent flf- in Fig. |I2|. For the bulk, mode-coupling theory 
predicts lim^^ = b |pD|] . Here, b is the von-Schweidler exponent which is related to 
1(f) by 70 = l/2a + 1/26 |22]-p4|. The exponents a (exponent of the "critical decay") and 
b are correlated with one another. A decrease or an increase of b entails a decrease or an 
increase of a. When applying these bulk predictions to the films the D-dependence of 0^ in 
Fig. [(3f(D = 20) < /3f(D = 5) < #f bulk ] suggests b(D = 20) < b(D = 5) < 6 bulk and thus 
70(.D = 2O) > 70 (.D = 5) > 70, bulk- While 7^(-D) > 70, bulk for all D, the results of table | do 
not confirm the expected order between the different film thicknesses. 

However, one should not conclude from the preceding discussion that a consistent ap- 
plication of MCT might not be possible. On the one hand, it is not clear whether the 
D-dependence of the finite-g result in Fig. [12| (/3^ at q — 6.9) is the same as in the limit 
q — > 00 (i.e., as for b). The corrections to the asymptotic behavior could depend on both q 
and D. On the other hand, the fitted value for 70 at D = 5 is the least certain because the 
lowest temperature simulated corresponds to a rather large T — T C (D = 5) ~ 0.045 contrary 
to T — T c ~ 0.025 for D = 20. Therefore, the simulations have to be extended to lower 
temperature to test whether there are clear deviations from the superposition of r q at low 
T in Fig. |TJ. The present quality of the fits would rather suggest that the 70-values for the 
different film thicknesses are very close to each other. 

Within the framework of MCT, quantities, such as the exponents a, b or 7 and the critical 
temperature, are determined by the thermodynamic properties of the glass former, in par- 
ticular by the static structure factor [[^HUl . The discussion of Fig. || already suggested that 
S(q) changes with decreasing D in a quite similar fashion as the bulk S(q) if the temperature 
is increased. Therefore, the question arises of whether it is possible to superimpose the bulk 
and film results for S(q) by comparing the data for the same T — T C (D). This would imply 
that the reduction of T c in the films is closely related to the fact that the development of 
the local packing, characteristic of the supercooled bulk, is shifted to lower temperature by 
the presence of the smooth walls. Figure [16] shows that such as superpostions for the same 
T — T C (D) is possible if the confinement is not too strong (here, for D > 10). The upper 
panel compares S(q) of the bulk with that of a film of thickness D — 10 for T — T c = 0.01. 
With the exception of the (slightly) different amplitude of the first maximum, both structure 
factors are identical over the whole g-range. The lower panel shows the same comparison 
for the bulk and films of thicknesses _D = 5, 10 and 20 at a larger distance from T c , i.e., for 
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T — T c = 0.05. While S(q) of the bulk and film still coincide for D > 10, this is no longer the 
case for the thinnest film studied (D — 5). The influence of the confinement on the packing 
structure of the system at this film thickness cannot be explained by a mere shift of the 
temperature axis. 

VI. SUMMARY 

Results of extensive MD simulations of thin (non-entangled) polymer films are presented, 
which focus on the influence of confinement on the sluggish dynamics of the system and in 
particular on the glass transition temperature. The film geometry is realized by introducing 
two perfectly smooth and purely repulsive walls. All simulations are carried out at constant 
normal pressure P^ext — 1- 

The static properties of the system show that chains close to the walls prefer a parallel 
alignment. However, when averaged over the whole film and all directions, the influence 
of the walls on the chain conformations becomes very weak. In particular, we find that 
the chains' radius of gyration, R 2 , does not depend much on film thickness. Even for the 
extreme case of D — 5, R 2 . lies only by 10% below the corresponding bulk value. 

On the level of the overall packing structure of the melt, we observe that the structure 
factor, S(q), of a film of thickness D resembles that of the bulk at a higher temperature. 
If the confinement is not too strong (D > 10), S(q), measured for a particular D at some 
T', almost coincides with the bulk result for that temperature T" which lies at the same 
distance to T c bulk as V to T C (D) [i.e. T" - T c bulk = V - T C {D)\. This indicates that T - T C (D) 
is a relevant parameter for our confined system. 

This static property of our model finds a counterpart in the dynamic behavior. Our main 
findings for the dynamics can be summarized as follows: (1) The relaxation of the super- 
cooled films is accelerated compared to the bulk so that characteristic temperatures, such as 
the mode-coupling critical temperature, T C {D), or the Vogel-Fulcher-Tammann temperature, 
T (D), decrease with decreasing film thickness. As we expect T < T g < T c , our results sug- 
gest that also T g (D) should decrease with decreasing film thickness. (2) The films exhibit 
several features predicted by mode-coupling theory, such as the space-time factorization 
property in the intermediate time window of the /3-process, time-temperature superposition 
of the a-relaxation, and a power-law increase of the a-relaxation time in a T-interval that is 
close, but not too close to T c . This implies that the implications of the cage effect, whose ap- 
proximate mathematical treatment leads to these predictions, also seems to be an important 
factor to understand the dynamics of the confined system in the supercooled state above T c . 
(3) A comparison of the mean-square displacements in direction parallel and perpendicular 
to the walls shows that not only the parallel motion, but also the perpendicular motion is 
accelerated compared to the bulk if the displacement is sufficiently smaller than the film 
thickness. However, for a given thickness, the parallel motion is always faster than that in 
transverse direction. In other words, the enhancement of the dynamics in more pronounced 
when relaxation processes parallel to the walls are considered. Due to the film geometry, it 
is clear that the long-time limit of the perpendicular mean-square displacements must be 
finite. We gave an expression which allows a computation of this limitting value from the 
density profile of the particle that corresponds to the MSD under consideration (i.e., inner 
monomer, center of mass, etc.). 
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TABLES 



D 


5 


7 


10 


15 


20 


bulk 


T 


0.204 ± 0.007 


0.253 ± 0.013 


0.288 ± 0.006 


0.297 ± 0.007 


0.308 ± 0.004 


0.328± 0.008 


T c 


0.305 ± 0.006 


0.365 ± 0.007 


0.390 ± 0.005 


0.405 ± 0.008 


0.415 ± 0.005 


0.450 ± 0.005 


7MSD 


2.5 ± 0.2 


2.4 ± 0.2 


2.1 ± 0.1 


2.2 ± 0.1 


2.1 ± 0.1 


1.84 ± 0.1 


70 


3.24± 0.08 


3.15 ± 0.1 


2.68 ± 0.08 


2.76 ± 0.1 


2.74 ± 0.1 


2.09 ± 0.07 



TABLE I. Survey of the VFT-temperature, Tq, the mode-coupling critical temperature, T c , 
and the critical exponents, 7msd an d 7^,, for different film thicknesses D and for the bulk. Tq 
was determined via fits to Eq. ( |l~7| ) both for the film and for the bulk. The critical temperature, 
T C (D), was obtained from fits to Eq. (14). In both cases, fits were done for the relaxation time 
extracted from mean-square displacements [see Eq. (|15|)]. T c bulk was known from previous analyses 
of incoherent and coherent scattering [ 17[ 4^j49]| . The same result for T c bulk is also obtained by 
applying Eq. ( |14| ) to the bulk-MSD's. The corresponding critical exponent 7msd is listed in the 



third row of the table. The critical exponent, 7^ (last row), is obtained from fits to Eq. (14), now 
using the a— relaxation times, r^.g, defined by (f> q (t=T g ) =0.3. In this case, T C (D) was kept fixed 
and only the prefactor and 7^ were varied. 
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FIGURES 




' 0.9 1.0 1.1 1.2 1.3 1.4 

r/G 

FIG. 1. Left panel: Lennard- Jones potential (LJ) versus the bond potential, the sum of the 
LJ and FENE potentials. The bond potential is shifted by -20 to lower values for the sake of 
comparison with the LJ potential. The minimum position of the bond-potential is smaller than 
that of the LJ-potential: The bond potential has its minimum at r ~ 0.96cr, whereas that of a 
pure LJ-potential lies at r = \pla. Due to the incompatibility of these length scales and due to the 
flexibility of our model (no bond angle or torsion potentials), one expects that cooling the system 
would not lead to crystallization, but maintain the amorphous structure typical of the liquid phase. 
This expectation is confirmed by the behavior of static structure factor [see lower panel of Fig. ||. 
Right panel: Schematic visualization of the local distortion of the structure due to the two intrinsic 
length scales of the model. 
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FIG. 2. Upper panel: A snapshot of the simulation box for a film of thickness D = 20 at 
T = 0.44. Note that the mode coupling critical temperature of the system at this film thickness 
is T C (D = 20) = 0.415 (see text and also |43|1 ). A temperature of T = 0.44 thus corresponds to the 
supercooled regime at this film thickness. To visualize the chain structure, a continuous color code 
is used, monomers of a given chain having the same color. Lower panel: Structure factor S(q) 
versus the modulus, q, of the wave vector for a film of thickness D = 10 at three temperatures: 
T = 0.42, 0.5, and T = l (the critical temperature of mode-coupling theory for this film thickness 
is T C (D = 10) ~ 0.39 ± 0.005, see Table ||). At low temperatures, the peak of the structure factor is 
more pronounced and its position is slightly shifted towards larger q. This is consistent with the 
fact that, at lower temperatures, the average interparticle distance decreases. Otherwise, S(q) does 
not develop sharp peaks for larger q at lower temperatures, thus indicating that the film remains 
amorphous. 
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FIG. 3. Upper panel: Chain structure factor, S c (q), for a film of thickness D = 20 averaged 
over the whole system. Two representative temperatures are shown: T = 1 (normal liquid state) 
and T = 0.44 (supercooled state) showing that S c {q) is insensitive to a temperature variation. 
S c (q) in the film is computed in the direction parallel to the walls, i.e. q = (q x ,Qy), Q = \q\, whereas 
q = (q x ,<ly,Qz) in the bulk. The solid line indicates the Debye function with Rg Debye = bulk 
[see Eq. (Q)]. The vertical dashed lines mark the q- values 2-ir/b (b = bond length) and 27r/i?g Ulk . 
The horizontal dashed line indicates the large g-limit of S c (q). Lower panel: Same as in the upper 
panel, but now evaluated for the bulk and for films of various thicknesses at T = 1 (normal liquid 
state). Apparently, S c (q) is hardly affected by the confinement. 
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FIG. 4. Upper panel: Components of the radius of gyration and of the end-to-end distance 
in directions parallel (R^ and R^ e and perpendicular (R^ j_ and R^ e j_) to the wall versus the 
distance, z, from the wall. R?Az) and R? Az) behave qualitatively similar. They develop a 

Sill ef =i I 

maximum close to the wall and then converge towards a constant (bulk) value in the film center 
(indicated by horizontal dashed lines). Both the bulk radius of gyration and the end-to-end distance 
are known from previous simulations: Rg bulk = 2.14, and i?g e bulk = 12.6 at T = 1 and pressure 
p=l ||16| . The parallel component is multiplied by a factor of 3/2 and the perpendicular one by 3. 
These factors take into account that compared to three independent directions in the bulk, there 
are only two independent coordinates in the parallel and one in the perpendicular directions in a 
planar system. Lower panel: Chain structure factor, S c (q), evaluated within layers of thickness 
of 1/4 monomer diameter for a film of thickness D = 5 and at the low temperature T = 0.35 
(supercooled state). Note that the mode-coupling critical temperature for this film thickness is 
T C (D = 5) = 0.305 db 0.005 [see Table [|. A layer close to the wall (•) and the layer at the film 
center (□) are compared to the average of S c (q) over the whole system (x). The chain structure 
in the film center coincides with the average behavior of the film. Thus, the deviations of the 
chain structure in the vicinity of the walls compared to the film center do not crucially influence 
the average over the whole system. The solid line gives the Debye function with Rg Debye = l.93 ; 
whereas the long dashed line corresponds to Rg Debye = 2.4 [see Eq. (^)J. The values of Rg Debye 
are the average radii of gyration calculated over all chains whose center of mass lies in the interval 
of thickness 1/4 around z = 1.125 and z = 2.375, respectively. The vertical dashed lines mark the 
(/-values 2ir/b (6 = bond length). The horizontal dashed line indicates the large-g behavior of S c (q) 
(i.e., S c (q -oo) = l). 
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FIG. 5. Upper panel: Parallel and perpendicular components of the radius of gyration (1.5 
R 2 I, and 3R 2 j_) and their sum, R 2 = R 2 ,, + R 2 ± , versus temperature. R 2 m and R 2 j_ have been 

o i 1 1 & ' " o ? 1 1 ^ ' ol 1 1 

computed as (time) averages over all chains in the system. Data are shown for two film thicknesses: 
D = 20 (open symbols) and D = 5 (the corresponding filled symbols). The factors 1.5 and 3 are 
introduced to simplify the comparison with the bulk value. They account for the different number 
of spatial components in the film for directions parallel (2) and perpendicular (1) to the walls 
compared to the bulk where there are 3 such components. Both in the case of D = 20 and of 
D = 5 chains prefer an alignment parallel to the walls. Due to the opposite changes of R 2 ,, and 
R 2 _|_ with respect to the bulk value, the effect of the confinement on the chain's radius of gyration, 
i?| = R 2 || + R 2 j_, is less strong and practically vanishes for D = 20. The horizontal dashed line 
indicates the bulk radius of gyration, R 2 bulk = 2.09 (at T = 0.5). Lower panel: Same quantities 
as in the upper panel, now versus film thickness, D, at a temperature of T = 0.5. (R^: connected 
circles, R^n- connected squares and Rg±- connected triangles). Again, the parallel orientation 
is favored by the walls. The overall radius of gyration, however, depends only slightly on film 
thickness. 
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FIG. 6. Upper panel: Structure factor, S(q), of the melt for the bulk and films of various 
thicknesses ranging from D = 5 (~ 3.5R g ) to D = 20 (~ 14i? g ) at a temperature of T = 1 (nor- 
mal liquid state). Note that the mode-coupling critical temperature of the system in the bulk is 
T c bulk = 0.45 [ 16 , 17 1 (see also Table |). S(q) in the film is computed by taking into account the 
coordinates parallel to the walls only, i.e. q = (q x ,Qy), Q=\q\- At T=l, the structure factor of the 
thickest film {D = 20) coincides with that of the bulk. For film thicknesses D < 10, the influence of 
the walls on S(q) is visible even at this high temperature. The peak value of S(q) decreases with 
stronger confinement (smaller D) and the position of the peak is slightly shifted towards smaller 
q, reflecting an increase of the (average) interparticle distance. Instead of decreasing the film 
thickness, these changes in S(q) can also be achieved by increasing the temperature only. So, the 
inverse film thickness plays a qualitatively similar role to that of the temperature. For the largest 
film thickness shown here, however, the structure factor of the film is identical to that of the bulk 
over a large g-range. Lower panel: The same as in the upper panel, now for a lower temperature 
of T = 0.46 (supercooled state). Contrary to T = 1, the structure factor of the film with D = 20 
no longer coincides with S(q) of the bulk. Thus, at lower temperatures, the confinement has a 
stronger impact on the packing structure of the melt. 
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FIG. 7. Comparison of the incoherent intermediate scattering function, (j) q {t), at T = 0.46 
( T buik = 0.45) for the bulk and films of various thicknesses ranging from D = 5 ~ 3.5R g to 
D = 20 ~ 14i? g . PN.cxt = 1 denotes the pressure of the simulation. The scattering functions 
are calculated at q = 6.9 (maximum of S(q), see Fig. §). The initial decay of <fi 8 q (t) is determined 
by free ballistic motion, i.e., (j) s q (t) = l — (ft s t) 2 /2 with Q, s = qyJk-QT . The dashed horizontal line at 
(f> q (t) = 0.1 and the vertical arrow on the curve for D = 5 indicate how the a-relaxation time, r q , is 
defined. It is the time value when the scattering function has decayed to 0.1. 



31 




32 



FIG. 8. Upper panel: Log-log plot of the mean-square displacement of innermost monomers, 
gi(t), at T = 1 (normal liquid state far above T^ n]k = 0.45). The figure compares the bulk data 
with the displacements measured parallel to the walls in films of different thicknesses ranging from 
D = 5 (« 3.5i? g ) to D = 20 (« 14i? g ). The film data were multiplied by 3/2 to account for the 
different number of spatial directions used to calculate g\ (t) (i.e., 3 directions for the bulk, but only 
2 for the films). The bulk end-to-end distance R 2 e (=12.3) and the average bond length b 2 = 0.922 
are indicated as horizontal dashed horizontal lines. The solid lines represent the behavior of gi(t) 
expected in different time regimes: ballistic at short times (gi(t) ~ t 2 ), diffusive at late times 
(gi (t) ~ t), and dominated by chain connectivity for times where g\(t) > b 2 (gi(t) ~ t x ; x = 0.63 = 
effective exponent in the bulk). For comparison, the MSD of chain's center of mass, gs (t), is 
also depicted: At short times, 53 (t) = Sk^Tt 2 /N p , where N p = 10 is the degree of polymerization 
(number of the monomers of a chain). Contrary to the innermost monomers, the motion of the 
chain's center of mass is not affected by chain connectivity. Therefore, 53 (t) does not show a 
t 63 -behavior, but continuously crosses over to the diffusive regime which is reached for 53 (t) > R 2 
(= 2.14 at T= 1). The effects of confinement are rather weak at this high temperature: gi(t) for 
D=10 and D = 20 are practically indistinguishable from that of the bulk. For strong confinement 
(D = 5 and D = 7), however, a dependence on film thickness is observed. Lower panel: Same as in 
the upper panel, now at a low temperature of T = 0.46 (supercooled state close to T c bulk = 0.45). 
The lowest horizontal line shows the plateau value 6r 2 c of a MCT-analysis (~ 0.054). Again, the 
solid lines represent the behavior of gi(t) expected in different time regimes: ballistic at short times 
(51(0 ~ t 2 )] diffusive at late times (g\(t) ~ t), and dominated by chain connectivity for times where 
<7i (t) > b 2 (gi(t) ~ t x ; x = 0.63 = effective exponent in the bulk). At this low temperatuure, the 
effect of the walls is much more pronounced compared to T = 1 (see the upper panel). In both 
panels PN,ext = 1 denotes the pressure of the simulation. 
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FIG. 9. Upper panel: Log-log plot of the mean-square displacement of all monomers, go(t), at 
T = 0.46 (supercooled state). The figure compares the bulk data with the displacements measured 
parallel (||) and perpendicular (_L) to the walls in films of different thicknesses: D = 5 (||: thick 
solid line, _L: thin solid line), D = 7 (||: thick dashed line, _L: thin dashed line) and D = 20 (||: thick 
dotted line, _L: thin dotted line). The film data were multiplied by 3/2 for the parallel direction 
and by 3 for the perpendicular one to account for the different number of spatial directions used to 
calculate go(t) (i- e -; 3 directions for the bulk, 2 for the parallel direction and 1 for the perpendicular 
direction in the film). The dashed horizontal line indicates the lavge-t limit of <7o,±(£) expected from 
Eq. @ for a film of thickness D = 5. Lower panel: Test of Eq. (||) for go,±(t = oo) and gi,_i_(t = oo) 
at T = l and T = 0.46 [go.j_(t ; T = 0.46) = connected circles, gi,±(t ; T = 0.46) = connected squares, 
<7o,±(* >T = 1) = connected diamonds and gi t ±(t ; T = 1) = connected triangles]. The film thickness 
is again D = 5. The density profiles of all monomers and of the innermost monomer are depicted 
in the inset. At the lower temperature of T = 0.46, both density profiles develop pronounced peaks 
close to the walls, whereas the density in the film center is rather decreased compared to T = 1. 
Larger transversal distances thus contribute more to the MSD in the z-direction. Therefore, we 
expect an increase of <7o,±(^ = °) an d gi,±(t = oo) at lower T [see also the text]. As seen in this 
panel, this expectation is indeed observed in the long time part of go,±(t) and 9i,j_(t), respectively. 
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FIG. 10. Test of the space-time factorization theorem [Eq. (|To|)1 for the incoherent intermediate 
scattering function <ffc(t) via the ratio Rq(t) defined by Eq. (|TTD . The three panels show the results 
for different film thicknesses and temperatures: D = 20, T = 0.46 (=T-T C = 0.045), 13 = 10, T = 0.42 
(=r-r c = 0.03), D = 5, T = 0.35 (=T-T C = 0.045). The smallest (q = l) and the largest (g=18.5) 
wave vectors are highlighted by a thick dashed line and a thick solid line, respectively. The thin 
solid lines in between refer to the following ^-values (from bottom to top): q = 1.8, 3, 4.3, 6.9, 7.1, 
9.5, 12.5, 16. The vertical dashed lines indicate the choices for t" (left line) and t' (right line). By 
definition, R s (t") = l and R s (t') = 0. 
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FIG. 11. Test of the time-temperature superposition principle [Eq. (|T^)1 for three differ- 
ent film thicknesses. Top panel: D = 20. Temperatures shown are from left to right: 
T = 0.46 (= T - T c = 0.045), 0.47, 0.48, 0.48, 0.5, 0.55 . Middle panel: D = 10. Temperatures 
shown are from left to right: T = 0.4 (=T-T C = 0.01), 0.42, 0.43, 0.44, 0.45, 0.46, 0.47, 0.5. Bottom 
panel: D = 5. Temperatures shown are from left to right: T = 0.35 (=T — T c = 0.045), 0.37, 0.38, 
0.39, 0.4, 0.41, 0.42. In all three panels the time axis is rescaled by the a-relaxation time, r q , 
defined by (p s (t = T q ) = 0.1, and the lowest temperature is depicted by a dashed line. 
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FIG. 12. Incoherent intermediate scattering function (j) q (t), calculated at the maximum of S(q) 
(q = 6.9), versus rescaled time t/r q for D = 5,20 and the bulk. r q is defined by (p s q (t = r q ) = 0.1. 
The temperatures shown roughly correspond to the same distance to T c , i.e., to T — T c ~ 0.045, 
for both the films and the bulk. In addition, the results of KWW-fits [Eq. (|l3|)] to the long-time 
behavior of (j) q (t) are depicted (D = 5: dashed line, D = 20: solid line, bulk: dot-dashed line). The 
corresponding stretching exponents are indicated in the figure. In all three cases, only data with 
t/T g > 0.05 were used for the KWW-Fit. 
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FIG. 13. Upper panel: Plot of t -1 / 7 versus T for D = 5. The relaxation time r was determined 
by Eq. (15) using the mean-square displacements (MSD) of inner-, end- and all monomers and of 
the chain's center of mass. The exponent 7 was first determined by fits to Eq. (|l~4|) with three open 
parameters. Since the results for all MSD's agreed within the error bars (7 = 2.5 ±0.2), 7 = 2.5 was 
used in the plot. The fits to Eq. (|lj|) are represented by straight lines. The intersection of these lines 
with the T-axis determines the critical temperature of the film (T C (D = 5) = 0.305 ± 0.005). Lower 
panel: Different representation of the same data. Now, r is plotted versus T—T c (T C (D = 5) = 0.305). 
The dashed line indicates the MCT-fit [Eq. (0)] using 7 = 2.5, whereas the solid line represents a 
fit to the VFT-equation [Eq. (PJJ)] with T (D = 5) = 0.204 ± 0.026. Both fits shown here were done 
for gi(t). 
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FIG. 14. Mode-coupling critical temperature, T C (D), and the VFT-temperature T${D) versus 
film thickness D. T C {D) was obtained from fits to Eq. (|14|). Similarly, the VFT-temperatures, 
Tq(D), and T bulk are results of fits to Eq. (|T7|). Note that, the glass transition temperature T g lies 
between these two temperatures: T c < T g < To. Thus, T g (D) lies somewhere between Tq(D) and 
T C (D). The solid line gives a suggestion for the possible form of T g (D). 
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FIG. 15. Upper panel: Relaxation times, r q , at the maximum of S(q) (q = 6.9) versus T — T c 
for various film thicknesses ranging from the smallest value studied (-D = 5~3.5i? g ) to the largest 
one (D = 20 « 14i? g ). T q was determined from the incoherent scattering function by requiring 
(jf q (t = Tq) = 0.3, whereas the critical temperatures, T C (D), were obtained from the analysis of the 
mean-square displacements, i.e. from power-law fits to r by Eq. (0). The lines represent fits of 
T q to Eq. ( |I4| ) in the regime where the data are linear: D = 5 (dash-dotted line), D = 10 (dashed 
line), bulk (solid line). The resulting values for 7<^(-D) are indicated in the figure. Upper panel: 
Using the critical exponents, j^D), obtained from the upper panel of this figure, r q l ^ 4 ' is plotted 
versus T — T C (D) for the film data. The lines are again fits to Eq. ([l4|). The deviations from the 
power law Eq. (0), visible for D = 10 at T — T C (D) = 0.01 in the upper panel, cannot be seen in 
this representation of the data. 
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FIG. 16. Upper panel: Collective structure factor, S(q), of the melt for a film of thick- 
ness D = 10 and for the bulk at temperatures T = 0.4 and T = 0.46, respectively. Note that 
T C {D = 10) = 0.39 and T^ nlk = 0.45. Thus, the chosen temperatures have the same distance from 
the corresponding critical temperature, i.e. 0A — T c (D = 10) = 0.46 — T c bulk = 0.01. Except the 
(slightly) different amplitude of the main peak, S(q) of the film and the bulk are essentially iden- 
tical when compared for the same difference to the critical temperature. Lower panel: Same as 
in the upper panel for a distance of 0.05 from T C (D), where data of many film thicknesses are 
compared to that of the bulk. Again, for film thicknesses D > 10, S(q) of the film overlaps with 
that of the bulk. For the extreme case of D = 5, however, deviations are not negligible. 
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